in 

q 

o 



A One- Vertex Decomposition Algorithm for Generating 
Algebraic Expressions of Square Rhomboids 

Mark Korenblit 1 and Vadim E. Levit 2 
1 Holon Institute of Technology, Israel 



(N 

o 

> 

O 1 

korenblit @hit .ac.il 
2 Ariel University Center of Samaria, Israel 



levit v@ariel .ac.il 



Abstract 



The paper investigates relationship between algebraic expressions and graphs. 
We consider a digraph called a square rhomboid that is an example of non-series- 
parallel graphs. Our intention is to simplify the expressions of square rhomboids and 
eventually find their shortest representations. With that end in view, we describe 
the new algorithm for generating square rhomboid expressions, which improves on 



\^ . our previous algorithms. 

1. Introduction 

A graph G = (V, E) consists of a vertex set V and an edge set E, where each edge corre- 
sponds to a pair (v, w) of vertices. A graph in which each edge has a weight associated 
with it is called a network. A graph G' = (V 7 , E') is a subgraph of G = (V, E) if V' Q V 
and E' C E. A graph G is a homeomorph of G' if G can be obtained by subdividing 
edges of G' with new vertices. We say that a graph G 2 = (V, E') is a square of a graph 
G = (V, E) if E' = {(u, w) :(u,w)eEV ((«, v) G E A (v, w) € E) for some v € V}. A 
two-terminal directed acyclic graph (st-dag) has only one source and only one sink. 

We consider a labeled graph which has labels attached to its edges. Each path 
between the source and the sink (a sequential path) in an st-dag can be presented by 
a product of all edge labels of the path. We define the sum of edge label products 
corresponding to all possible sequential paths of an st-dag G as the canonical expression 
of G. An algebraic expression is called an st-dag expression (a factoring of an st-dag in 
[4]) if it is algebraically equivalent to the canonical expression of an st-dag. An st-dag 
expression consists of literals (edge labels), and the operators + (disjoint union) and • 
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Figure 1.1: A series-parallel graph. 

(concatenation, also denoted by juxtaposition). An expression of an st-dag G will be 
hereafter denoted by Ex(G). 

We define the total number of literals in an algebraic expression as the complexity 
of the algebraic expression. An equivalent expression with the minimum complexity is 
called an optimal representation of the algebraic expression. 

A series-parallel graph is defined recursively so that a single edge is a series-parallel 
graph and a graph obtained by a parallel or a series composition of series-parallel graphs 
is series-parallel. As shown in [4] and [12], a series-parallel graph expression has a rep- 
resentation in which each literal appears only once. This representation is an optimal 
representation of the series-parallel graph expression. For example, the canonical expres- 
sion of the series-parallel graph presented in Figure 1.1 is abd+abe + acd+ace + fe + fd. 
Since it is a series-parallel graph, the expression can be reduced to (a(b + c) + /) (eZ + e) , 
where each literal appears once. 

A Fibonacci graph [10] has vertices {1, 2, 3, . . . , n} and edges {(v, v + 1) | v = 1, 2, . . . , 
n — 1} U {(v, v + 2) | v = 1, 2, . . . , n — 2}. As shown in [5], an st-dag is series-parallel if 
and only if it does not contain a subgraph which is a homeomorph of the forbidden sub- 
graph positioned between vertices 1 and 4 of the Fibonacci graph illustrated in Figure 
1.2. Thus, a Fibonacci graph gives a generic example of non-series-parallel graphs. 
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Figure 1.2: A Fibonacci graph. 

Mutual relations between graphs and expressions are discussed in a number of 
works. Specifically, [18], [19], and [25] consider the correspondence between series- 
parallel graphs and read-once functions. A Boolean function is defined as read-once if 
it may be computed by some formula in which no variable occurs more than once (read- 
once formula). On the other hand, a series-parallel graph expression can be reduced to 
the representation in which each literal appears only once. Hence, such a representation 
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of a series-parallel graph expression can be considered as a read-once formula (boolean 
operations are replaced by arithmetic ones). 

Problems related to computations on graphs have applications in different areas. 

Specifically, many network problems, which are either intractable or have compli- 
cated solutions in the general case are solvable for series-parallel graphs. For example, 
some efficient algorithms for flow problems on series-parallel networks are presented in 
[2], [3], [11], [26]. Papers [1], [6], [16], [17] consider sequencing and scheduling in relation 
to precedence series-parallel constraints. Linear algorithms for reliability problems on 
series-parallel networks are presented in [24], [27]. 

An expression of a homeomorph of the forbidden subgraph belonging to any non- 
series-parallel st-dag has no representation in which each literal appears once. For 
example, consider the subgraph positioned between vertices 1 and 4 of the Fibonacci 
graph shown in Figure 1.2. Possible optimal representations of its expression are 
o-i (02O3 + ^2) + ^103 or (a±a2 + &i) 03 + ai&2- For this reason, an expression of a non- 
series-parallel st-dag can not be represented as a read-once formula. However, for ar- 
bitrary functions, which are not read-once, generating the optimum factored form is 
NP-complete [28]. 

The problem of factoring boolean functions into shorter, more compact formulae is 
one of the basic operations in algorithmic logic synthesis. In logic synthesis, one standard 
measure of the complexity of a logic circuit is the number of literals. Computation time 
also depends on the number of literals. Some algorithms developed in order to obtain 
good factored forms are described in [8], [9]. 

A symbolic approach to scheduling of a robotic line is considered in [15]. The method 
uses the max-algebra tools and allows the shortest-path problem to be interpreted as the 
computation of the st-dag expression. The complexity of this problem is determined by 
the complexity of the st-dag expression. For a robotic line represented by a Fibonacci 
graph, the proposed algorithm generates the processing sequence in polynomial time. 

A method for automated composition of algebraic expressions in complex business 
process modeling based on acyclic directed graph reductions is introduced in [20]. The 
method transforms business step dependencies described by users into digraphs and 
finally generates algebraic expressions. If a graph is not series-parallel, the algorithm 
checks potential structural conflicts whose presence complicates certain aspects, such 
as execution control and system scalability. In this case, the expression generation may 
require exponential time. 

A lot of problems have polynomial algorithms on planar graphs, i.e., the graphs that 
can be drawn in the plane without any edges crossing. For example, some polynomial- 
time algorithms for the exact computation of network reliability applied to planar graphs 
are surveyed in [21]. In [22] a polynomial algorithm for the connectivity augmentation 
problem is presented. This problem has applications in upgrading telecommunication 
networks to be invulnerable to link or node failures. There are also many practical 
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applications with a graph structure in which crossing edges are a nuisance, including 
design problems for circuits, subways, utility lines [7]. 

In our work we consider expressions with a minimum (or, at least, a polynomial) 
complexity as a key to generating efficient algorithms on distributed systems. 

In [12] we presented an algorithm, which generates the expression of O (n 2 ) com- 
plexity for an n-vertex Fibonacci graph. 

In this paper we investigate a non-series-parallel st-dag called a square rhomboid 
(Figure 1.3). This graph looks like a planar approximation of the square of a rhomboid, 
which is a series composition of rhomb graphs. A square rhomboid consists of the same 
vertices as the corresponding rhomboid. However, edges labeled by letters a, b, and c 
(see Figure 1.3) are absent in a rhomboid. Geometrically, a square rhomboid can be 
considered to be a "gluing" of two Fibonacci graphs, i.e., it is the next harder one in a 
sequence of increasingly non-series-parallel graphs. 




N-l 



Figure 1.3: A square rhomboid of size n. 
The set of vertices of an N- vertex square rhomboid consists of middle (basic), 



— g— upper, any* — g- 

denoted in formulae by x and x, respectively. The square rhomboid (SR for brevity) 
including n basic vertices will be denoted by SR(n) and will be called an SR of size n. 

Our intention in this paper is to simplify the expressions of square rhomboids and 
eventually find their optimal representations. With that end in view, we present an 
algorithm based on a decomposition method. 



and 



N-l 



lower vertices. Upper and lower vertices numbered x will be 



2. A One- Vertex Decomposition Method (1-VDM) 

The method is based on revealing subgraphs in the initial graph. The resulting expres- 
sion is produced by a special composition of subexpressions describing these subgraphs. 

For a non-trivial SR subgraph with a source p and a sink q we choose any decompo- 
sition vertex i (p < i < q) located in the basic group of a subgraph. We conditionally 
split each SR through its decomposition vertex (see the example in Figure 2.1). 
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Figure 2.1: Decomposition of a square rhomboid of size 7 at vertex 4. 



Two kinds of subgraphs are revealed in the graph in the course of decomposition. 
The first of them is an SR with a fewer number of vertices than the initial SR. The 
second one is an SR supplemented by two additional edges at one of four sides. Possible 
varieties of this st-dag (we call it a single-leaf square rhomboid and denote by SR) which 
are subgraphs revealed from an SR in Figure 2.1, are illustrated in Figure 2.2(a, b, c, 
d). Let SR(n) (an SR of size n) denote an SR including n basic vertices. 

We denote by E(p, q) a subexpression related to an SR subgraph with a source p 
and a sink q. We denote by E(p,q), E(p,q), E(p,q), E(p,q) subexpressions related to 

SR subgraphs of Figure 2.2(a, b, c, d, respectively) with a source p and a sink q, a 
source p and a sink q, a source p and a sink q, and a source p and a sink q. 




Figure 2.2: Varieties of a single- leaf square rhomboid of size 3. 

Any path from vertex 1 to vertex 7 in Fig. 2.1 passes through decomposition vertex 
4 or through edge C3 or through edge 03. Therefore, E(p, q) is generated by the following 
recursive procedure (decomposition procedure): 
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1. case q = p : E(p, q) <— 1 

2. case q = p + 1 : E(p, q) ^b p + e 2 p-ie 2 p + d 2 p-id 2 p 

3. case g > p + 1 : choice(p, q, i) 

4. £(p, g) <- E(p, i)E(i, q) + E(p, T^l^EfJ, q)+ 

E(p, i - l )g i _i.E(i,g) 

The expression related to a one- vertex .Si? is denned formally as 1 (line 1). Line 2 de- 
scribes the expression of SR(2) which is a trivial subgraph. The procedure choice(p, q, i) 
in line 3 chooses an arbitrary number i on the interval (p, q) so that p < i < q. This 
number is a number of a decomposition vertex in a basic group of the SR subgraph 
positioned between vertices p and q. A current subgraph is decomposed into six new 
subgraphs in line 4. Subgraphs described by subexpressions E(p, i) and E(i, q) include 
all paths from vertex p to vertex q passing through vertex i. Subgraphs described by 
subexpressions E(p, i — 1) and E(i, q) include all paths from vertex p to vertex q passing 
via edge Cj_i. Subgraphs described by subexpressions E(p, i — 1 ) and E(i,q) include all 
paths from vertex p to vertex q passing via edge ctj-i. 

E(l,n) is the expression of the initial SR of size n. Hence, the decomposition 
procedure is initially invoked by substituting 1 and n instead of p and q, respectively. 

We now consider the algorithm that generates the algebraic expression for a square 
rhomboid using 1-VDM as its base. 

3. A One- Vertex Decomposition Algorithm (1-VDA) 

The algorithm is based on decomposition of the initial graph and all kinds of subgraphs 
to new subgraphs. 

We conjecture that the representation with the minimum complexity of Ex(SR) 
derived by 1-VDM is generated when the decomposition vertex is a middle vertex in the 
basic group of the SR. That is, the number i of the decomposition vertex for a current 
SR subgraph which is positioned between vertices p and q is chosen as ( j" 2 ^ 2 ] or 
L 2 ?] )■ I n °t ner words, i is equal to for odd q — p + 1 and i equals q+ 2~ 1 or q+ P +1 
for even q — p + 1 . 

An SR subgraph is decomposed through a decomposition vertex selected in its basic 
group into six new subgraphs in the same way as an SR (see the examples in Figure 
3.1). The decomposition vertex is chosen so that the location of the split is in the middle 
of the subgraph. 

Three kinds of subgraphs are revealed in an SR in the course of decomposition. The 
first and the second of them are an SR and an SR, respectively. The third one is an SR 
supplemented by two additional pairs of edges (one pair is on the left and another one 
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Figure 3.1: Decomposition of single- leaf square rhomboids by 1-VDA. 



is on the right). Possible varieties of this st-dag (we call it a dipterous square rhomboid 

and denote it by SR) are illustrated in Figure 3.2(a, b, c, d). We define subgraphs 

illustrated in Figure 3.2(a, c) as trapezoidal SR graphs and subgraphs illustrated in 

Figure 3.2(b, d) as parallelogram SR graphs. Let SR(n) (an SR of size n) denote an 

SR including n basic vertices. 

We denote by E(j>,q), E(p,q), E(p,q), E(p,q) subexpressions related to subgraphs 
of Figure 3.2(a, b, c, d, respectively) with a source p and a sink q, a source p and a sink 
q, a source p and a sink q, and a source p and a sink q. 




Figure 3.2: Varieties of a dipterous square rhomboid of size 2. 



An SR subgraph is decomposed into six new subgraphs in the same way as an SR 
and an SR (see the examples in Figure 3.3). The number i of the decomposition vertex 

in the basic group for a current SR subgraph which is positioned between vertices p 
and q, is chosen as 9+ ?? +1 ( 



2 



or 



q+p+l 
2 



In the course of decomposition, two kinds of subgraphs are revealed in an SR. 

They are an SR and an SR. Therefore, no new kinds of subgraphs are formed and 
subexpressions related to all subgraphs can be generated. 

Hence, Ex(SR) is computed by the following recursive relations: 
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Figure 3.3: Decomposition of dipterous square rhomboids by 1-VDA. 



1. E(p,p) = 1 

2. E(p, p) = e 2p -i 3. E(p,p) = d 2p -i 4. E{p,p + l) = e 2p 5. E(p,p + 1) = d 2p 
6. E(p, p+ 1 ) = c p + e 2p e 2p+1 7. E(p, p + 1 ) = e 2p (i2p+i 

8. E(p,p+ 1) = d 2 pe2 P +i 9. E(p, p+1 ) = a p + d 2p d 2p+ i 

10. E(p, p+ 1 ) = 6 P + e 2 p_ie 2 p + d 2p -\d 2p 

11. £(p,p+l) = (ftp + d 2p -id 2p )e 2p+1 + e 2p _i(c p + e 2p e 2p+ i) 

12. E(p, p + 1 ) = (6p + e 2 p_ie 2 p)(i 2 p + i + d 2 p_i(a p + d 2p d 2p+1 ) 

13. E(p,p + 2) = (cp + e 2 pe 2 p+i)e 2p+2 + e 2p (b p+ i + (^+1^+2) 

14. E(p, p + 2 ) = (dp + d 2p d 2p+ i)d 2p+2 + d 2p (b p+1 + e 2 p + ie 2 p +2 ) 

15. E(p,p + 2) = e 2p (b p+ i + d 2p+ id 2p+2 )e 2p+3 + (a p + d 2p d 2p+ i)(a p+ i + ^+2^+3) 

16. E(p, p + 2 ) = e 2 p(6p + id 2 p +3 +(i 2 p + i(ap + i+(i 2 p +2 d 2 p + 3))+(cp+e 2 pe 2 p + i)e 2 p +2 d 2 p + 3 

17. E(p,p + 2) = d 2p (bp +1 e 2 p +3 +e 2p+ i(c p+ i+e 2 p +2 e 2 p +3 ))+(a p +d 2p d 2p+ i)d 2 p +2 e 2 p +3 

18. E(p, p + 2 ) = d 2p (b p+ i + e 2p+ ie 2p+2 )d 2p+3 + (c p + e2pe 2 p+i)(cp + i + e 2p+2 e 2p+3 ) 

19. £(p, 5) = £7(p, g)+£(p, t^l)c i _i J E7(» ) i - lK-i£(j, g), i = [2+E] 

(«Z>P+1) 

20. £(p,g) = E(p,i)E(i,q)+E(p,i- l)c i - 1 E(i,q)+E(p,i^ 1 l)a i - 1 E(i,q), i = 

(q>P+l) 

21. £(p, g) = £7(p, 9)+^(p, i - l)Q-i£(*, ?)+S(p, i - l)aj-iE(i, q),i= [2+£] 
(9>P + 2) 

22. £(p, g ) = E(p,i)E(i,q)+E(p,i- l)c i - 1 E(i,q)+E(p,i^ 1 l)a i - 1 E(i,q), i = [2±2] 
(9>P+2) 

23. £(p,g) = E(p,i)E{i,q)+E{p,i- 1)0^(1, q)+E(p,i^ 1 l)a i - 1 E(i,q), i = 2±|±1 
(9>P+2) 

24. £(p,g) = E(p,i)E(i,q)+E(p,i- l)c i - 1 E(i,q)+E(p,i^ai- 1 E(i,q), i = 2±f±i 
(9>P+2) 

25. £(p,g) = E(p,i)E(i,q)+E{p,i - l)c i - 1 E(i,q)+E(2,i=r)a i - 1 E(i,q), i = 
(9>P+2) 

26. £(p,g) = E(p,i)E(i,q)+E(p,i- 1)0^(1, q)+E{p,i^ 1 l)a i - 1 E(i,q), i = 2±|±i 

(q>P + 2) 

27. £(p,g) =E(p,i)E(i,q)+E(p,i-l)c i - 1 E(i,q)+E(p,i^a i -iE(i,q), i = 2±2 
(<7>p+l). 



Lemma 3.1. Complexities of expressions Ex ^trapezoidal SR(n) ) and 
Ex I parallelogram SR(n) I derived by 1-VDA are equal for n > 2. 



Subgraphs of sizes 1 and 2 are trivial and their expressions are in lines (1 - 18). 

Specifically, expressions of subgraphs SR(2) and SR(2) are presented in the minimum 
factored form (lines 11 - 18). 

For example, the expression of the square rhomboid of size 3 derived by 1-VDA is 

(bi + eie 2 + + e3e 4 + d 3 d 4 ) + e\C\e± + d\a\d<±. 

It contains 16 literals. 

The following lemma results from lines (23 - 26) of 1-VDA. 



The following proposition results from Lemma 3.1 and from lines (1 - 27) of 1-VDA. 

Proposition 3.2. The total number of literals T(n) in the expression Ex(SR(n)) de- 
rived by 1-VDA is defined recursively as follows: 

1) T(l) = 0; 2) f (1) = 1; 3)T pr (l) = 2; 4) f Jr (l) = 3 

5) T{2) = 5; 6)f (2) = 8; 7) f pr {2) = 12; 8) f tr (2) = 11 

9) f (3) = 22; 10) f(3) = 28; 11) f (4) = 47; 12) f (4) = 60 

13) f (5) = 79; 14) f (5) = 92; 15) f (6) = 132; 16) f(6) = 50 

17) T(n) = T T (Lf J +1)+2T ([2] - l) +2? ([fj) +2 (n > 2) 

18; T(n) = T ( [f J + 1) + f ( [J] ) + 2T ( [f J ) + 2T ( [| ] - l) + 2 (n > 6) 

19; T(n) = f ( 1 ) + f ( Lf J + 1) + 2T ( [5 1 - 1) + 2T ( [f J ) + 2 (n > 6), 

wiere T(n), T pr (l), T tr (l), and T(n) are the total numbers of literals in Ex(SR(n)), 



Ex ^parallelogram SR (1) J , Ex ^trapezoidal SR (1) J , and Ex ySR (n) J , respectively. 

Hence, the expression Ex(SR(n)) derived by 1-VDA consists of six subexpressions 
related to six revealed subgraphs of size n' ~ f and two additional literals. Each 
of these subexpressions is constructed in its turn from six subexpressions related to 
six subgraphs of size n" ~ j and two additional literals, etc. Thus, by the master 
theorem, T(n) = 0(n log 2 6 ), i.e., 1-VDA provides the representation of Ex{SR{n)) 
with a polynomial complexity. 

It is of interest to obtain exact formulae describing complexity of the expression 
Ex(SR(n)) derived by 1-VDA. We attempt to do it for n that is a power of two, i.e., 
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n = 2 k for some positive integer k > 2. Statements (17 - 19) of Proposition 3.2 are 
presented for n = 2 k (k > 3) as 



' T (n) = T (§ + 1) + T (§) + 2T (§) + 2T (§ - l) + 2 

f (n) = T (f + l) + 3T (§) + 2? (f - 1) + 2 
, f (n) = f + 1) + f (§ ) + 2T (§) + 2? - 1) + 2 , 



(3.1) 



respectively. 

After a number of transformations, the following explicit formulae for simultaneous 
recurrences (3.1) are obtained by the method for linear recurrence relations solving [23]: 



1 54 

T( n ) = — n 10 ^ 6 
V ' 135 

(Z, \ 154 

T(B) = 135" 

^/ \ 154 
T(n) = —n 



1 log 2 3 _ 2 

27 5 



2 « _|_ 

logo 6 



19 



log 2 3 _ 



27 

a6 + -:-' 



(3.2) 



l°g2 3 

27" 5' 



4. Comparison of 1-VDA with Other Algorithms 

Another decomposition method (called 2-VDM) for generating algebraic expressions 
of square rhomboids is presented in [13] and [14]. The method consists in splitting 
a square rhomboid through two decomposition vertices one of which belongs to the 
upper group and another one to the lower group. The following algorithms based 
on this method are considered in [13] and [14]: full decomposition algorithm (FDA), 
combined decomposition algorithm (CDA) and improved FDA (IFDA). Complexities of 
generated by these algorithms expressions are also O (n log26 ). The values of T(n) for 
these algorithms and for 1-VDA are presented in Table 4.1 (4 < n < 10). 



n 


T(n), FDA 


T(n), CDA 


T(n), IFDA 


T(n), 1-VDA 


4 


47 


43 


43 


41 


5 


110 


102 


100 


66 


6 


173 


161 


157 


119 


7 


252 


236 


228 


172 


8 


331 


311 


299 


247 


9 


520 


488 


470 


322 


10 


709 


665 


641 


439 



Table 4.1: Complexities for FDA, CDA, IFDA, and 1-VDA. 



One can see that 1-VDA is more efficient than FDA, CDA, and IFDA. 

10 



The numerical advantage of 1-VDA in contrast to other algorithms follows also from 
presented in [13] and [14] explicit formulae analogous to (3. 2). The coefficient of the 
leading term of the formulae - n log 2 6 is equal to || « 1.76 for FDA, y|| 1.68 for 
CDA, fi| « 1.57 for IFDA, and ±|| « 1.14 for 1-VDA. 

' 135 ' 135 

5. Conclusions 

The decomposition method divides all subgraphs, which are revealed in the current 
recursive step into a constant number of subgraphs of proportionally decreasing sizes. 
The existence of a decomposition method for an st-dag G is a sufficient condition for the 
existence of a representation of Ex (G) with a polynomial complexity. The complexity 
depends on the number of revealed subgraphs in each recursive step. 
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